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In this paper we report results from a systematic strong-coupling expansion of the Bose-Hubbard 
model in one and two spatial dimensions. We obtain numerically exact results for the structure 
factor and the spectrum of single particle and single hole excitations in the Mott insulator. This 
enables the determination of the zero-temperature phase diagram and the location of the critical 
endpoints of the Mott lobes. In one dimension we confirm the occurrence of reentrance behavior 
from the compressible to the insulating phase in a region close to the critical point. 

PACS numbers: 05.30. Jp, 05.70.Jk, 67.40.Db 



I. INTRODUCTION 

Quantum phase transitions in strongly correlated sys- 
tems have attracted a lot of interest in recent years. Here 
we focus on systems of interacting bosons. Physical real- 
isations include Josephson junction arrays, granular and 
short-correlation-length superconductors, flux-lattices in 
type-II superconductors and possibly in the future ultra- 
cold atoms in a periodic potential. 

The minimal model containing the key feature of com- 
petition between kinetic and potential energy is the Bose- 
Hubbard hamiltonan. It is given by 



H 



-t 



U 



where the h\ and hi are bosonic creation and annihilation 
operators, fii = h\hi is the number of particles on site i, t 
the hopping matrix element, [/ > the on-site repulsion 
and /i the chemical potential. 

This model has been ^W^4p|4| flfjEviously by Quantum 
Monte Carlo simulationsElu'tra'Bflu in one and two spa- 
tial dimensions. Recently the one-dimensional case was 
also investigated using the density-matrix renormalisa- 
tion group (DMRG). This study found indications for 
reentrence from the superfluid to the Mott-insulator for 
certain parameter values. 

The zero-temperature phase diagram can be under- 
stood by starting from the strong-coupling or "atomic" 
limit where the the kinetic energy term vanishes [t = {)). 
In this limit the ground state is an insulator with a fixed 
number no of particles per site and the wavefunction is 
given by 



2 ^ 



I) 



(1) 



FO^Mott 



N 

n 



1 



|0> 



(2) 



^Mott/^ = U^noino - 1) - fino 



(3) 



Single charge excitations in the atomic limit are created 
by adding or removing a particle onto or from a particular 
site i: 
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The energy of these states relative to the ground state is 
given by 
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for particle and hole excitations respectively. For certain 
critical values of the chemical potential fi^^ = Uhq this 
energy vanishes and the system becomes compressible. 
If the hopping matrix element t is finite the range of the 
chemical potential for which the system is incompress- 
ible decreases. At some critical value tc of the hopping 
matrix element the Mott insulator will completely disap- 
pear. The resulting phase diagram shows Mott insulating 
regions, usually referred to as Mott lobes, surrounded by 
a compressible phase. 

This phase is characterized by the following real space 
correlation function: 



S, = {b]bo + b,bl) 
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In spatial dimension D > 1 the correlations are long 
ranged and the system is a superfluid, while in one di- 
mension it is a Luttinger Liquid with Sj decaying alge- 
braically for large distances. In this article we will inves- 
tigate the Fourier transform 
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and it's second moment 
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Assuming a lorentzian shape of >S'(q) near the ordering 
wave vector q = the correlation length, ^, is given by 
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We used strong coupling expansions in t/U to calcu- 
lated the ground state energy, Eq, the spectrum of single 
charge excitations Ep^j-f^^i^/f^gi^ (q) , the correlation func- 
tion ^(q) and the correlation length ^. The resulting 
series were analysis using Pade extrapolation techniques. 

In the next section we give a technical outline of how to 
calculated the ground state energy, Eq^ the spectrum of 
single charge excitations EparUcie/hoieil), the correlation 
function S{q) and the correlation length ^ by strong cou- 
pling expansions in t/U. We use methods that had been 
developed over the last decade but so far were mainly 
applied to spin systems. 

Section III focusses on the ground state and the ex- 
citated states. We present the dispersion of the single 
particle excitations and the dependence of the particle- 
hole gap on the hopping amplitude t. It is found that 
the series converges rather fast provided t is smaller then 
the critical value tc- In this regime, the Mott insulator, 
we obtain the numerical exact single charge excitation 
spectrum. 

In the following section the phase diagram in two di- 
mensions is discussed. Again due to the rapid conver- 
gence of the series it is possible to determine numerically 
exact the boundary of the Mott lobes. Only in a small 
region close to the critical point tc does the convergence 
break down and series extrapolation techniques have to 
be applied. These also work extremely well and it is thus 
possible to determine the critical point with high accu- 
racy. We also determine the critical exponents and find 
that they agree with the field theoretic predictions that 
this D-dimensional quantum system is in the universality 
class of the the D + 1-dimensional classical xy-model. 

Section V deals with the special case of a one- 
dimensional system. As already mentioned the compress- 
ible phase exhibits only quasi longrange order. Also the 
phase transitions at the top of the Mott lobes changes 
qualitatively and becomes of the Bereshinski-Kosterlitz- 
Thouless (BKT) type. Our analysis of the phase diagram 
finds a reentrance behaviour, i.e. for a certain range of 
values of the chemical potential the system has not just 
one transition from the Mott state to the compressible 
phase but upon further increase of the hopping strength t 
becomes insulating again until a second transition to the 
liquid occurs. This scenario agrees with recent DMRG 
calculations. We give a simple intuitive explanation for 
this surprising observation. 

Conclusion 



II. SERIES EXPANSION TECHNIQUE FOR 
EXCITED STATES 



Strong coupling perturbation theoryi-can be formulated 
as a linked cluster expansion, see e.gJ3. For the ground 
state energy this is basically nondegenerate Rayleigh- 
Schrodinger perturbaton theory performed for every clus- 
ter contributing up to a certain order in the expansion. 
In order to obtain correlation function one has to formaly 
add source terms. 
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to the hamiltonian and expand the energy to linear or- 
der in the J^j-. Differentiating with respect to the source 
fields will give the strong coupling expansion for the cor- 
relation functions. This procedure is necessary, because 
it is not possible to evaluate the wave function directly 
by a linked cluster expansion. This is due to the fact 
that the linked cluster theorem applies only to physical 
observables that are additive when the system seperates 
into disconnected part. 

A systematic strong coupling expansion of the energy 
of the charge excitations complicated due to the high 
degeneracy. The problem how to write down a linked 
cluster expansion fog degenerate states was solved only 
recently by GelfandB. The idea is to construct pertur- 
batively an effective Hamiltonian iJfJ in the subspace of 

the degenerate states |no; *)part/hoio ^ similarity trans- 
formation 



with 
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where Greek indices run over states in the full Hilbert 
space while Latin indices are restricted to the degenerate 
manifold of single particle and single hole states (^ and 
(H) respectively. Then the linked cluster theorem applies 
to Hfj{t) — Euottit)- In the case of a homogeneous sys- 
tem depends only on the difference of indices i — j 
and is easily diagonalised by a Fourier transform. This 
way one can determine the full dispersion E{q;, t, ji) of 
the charge excitations. In many ways the linked cluster 
expansion is similar to a exact diagonalization study of 
small systems - however in the linked cluster expansion 
it is possible to remove all finite size effects in each order 
and one obtains the full dispersion in the thermodynamic 
limit. 

We close this section with some technical remarks. Ex- 
pansions were performed up to 13*'' order for the tri- 
angular lattice with 12,253 clusters contributing. The 
Hilbert-space of the largest cluster had 37,442,160 states 
in the the single particle sector with filling no = 1. Up 
to 13*'' order in perturbation theory just a fraction of 
these states (788,238) actually contribute to the expan- 
sion. It was only after we restricted the calculation to this 
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FIG. 1. Dispersion of the single-particle and single-hole 
excitations of the square lattice Bose-Hubbard model at 
t/U = 0.055 



FIG. 2. Correlation function S{k) of the square lattice 
Bose-Hubbard model at t/U = 0.055 



subset of relevant states that it was possible to perform 
this work. Any significant improvement is unlikely be- 
cause the Hilbert-space and the number of clusters both 
increase exponentially with the order of the expansion. 

III. GROUND STATE AND SINGLE CHARGE 
EXCITATIONS 

We will first discuss the two dimensional case. We 
investigated both the square and triangular lattice and 
calculated the series for occupation numbers uq — 1 and 
no = 2 up to IS*'' and 10*'^ order respectively. 

The spectrum i?(q; t, fi) takes on the form 

^^part(q;i,Ai) epart(q;i) - A* (14) 
£^hoio(q; t, = -ehoic(q; t) + ^ (15) 

in complete analogy to eqns. and (|7|). For posi- 
tive values of the hopping matrix element t the small- 
est (largest) eigenvalue in the particle (hole) sector is al- 
ways located at wavevector q = 0. The upper and lower 
phase boundaries of the Mott phase are thus given by 
A<uppor(i) = epart(q = 0; t) and /xiowor(0 = ehoio(q = 0;t), 
respectively. 

As a consequence the single charge gap A(t) = 
epart(q = 0; t) — ehoio(q = 0; t), determines also the width 
A*upper(i) — /Slower (i) of the insulatlng region. With in- 
creasing hopping t the distance between the upper and 
lower boundary decreases until finally at some critical 
value, tc, the energy to remove a particle and the en- 
ergy to add a particle become degenerate and the Mott 
insulator vanishes altogether. 



The dispersion of the particle and hole excitations for 
?io = 1 on the square lattice and triangular lattice is 
shown in Figs. |^ and ||. The different shape of the two 
curves refiects the particle-hole asymmetrie of the model 
hamiltonian (|^). The series were found to converge very 
fast. Both Figures were obtained by summation of the 
13*'^ order series. They turned out to be almost indistin- 
guishable from the result of the 10-term series even for 
t/U = 0.055 (square lattice) and t/U = 0.035 (triangular 
lattice) which is not far from the critical endpoints t^ of 
the first Mott lobe for these two lattices. The particle 
and hole excitations both have a pronounced extremum 
at wavevector q = and are separated by a gap A. For 
values of the chemical potential ^ in this range all single 
charge excitations are gapped and the system is insulat- 
ing. Coefficients of the series expansion of A are given in 
Tab.(l). 

Figs.^ and|^ show the correlation function S{q) for the 
same sets of parameters as for the dispersion. Although 
the system is still in the Mott phase there is already a 
prounced peak at wavevector q = 0. 

IV. 2D PHASE DIAGRAM 

The phase diagram can be obtained from the series for 
the single charge exciations. Outside the critical region 
the values of the chemical potential for which the system 
becomes compressible can be directly calculated from the 
series. Only near the critical point is it necessary to use 
extrapolatation methods. We applied a Fade analysis 
and obtained the phase diagrams shown in Figs. I and I 
for the square and triangular lattices. With the excep- 
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FIG. 3. Dispersion of the single-particle and single-hole 
excitations of the triangular lattice Bose-Hubbard model at 
t/U = 0.035 
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FIG. 4. Correlation function S{k) of the triangular lattice 
Bose-Hubbard model aX t/U — 0.035 



tion of the lowest approximand all others turn out to be 
almost indistinguishable again indicating rapid conver- 
gence. 

The bare series shows a remarkable fast convergence 
for not too large values of the expansion parameter. In 
each order k of the expansion A vanishes at some effective 
critical value tc{k) with a corresponding effective fJ.c{k). 
Plotting tc{k) and fidk) vs. 1/fc one finds again a rapid 
convergence as shown in Figs.(|^) and (^. Extrapolation 
to n — > CX3 allows to determine accurately the critical 
point: 

tc = 0.05974 ± 0.00004 (square lattice) (16) 
tc = 0.037785 ± 0.000005 (triangular lattice) (17) 

The accuracy for the triangular lattice is almost an order 
of magnitude better then for the square lattice. This is 
because the series for the triangular lattice has a mono- 
tonic convergence while the square lattice shows an oscil- 
latory behaviour between even and odd terms as can be 
seen when comparing Figs.(0) and (|). The convergence 
of the data is significantly accelerated when considering 
the difference between neighbouring data points. Fitting 
each pair of points tc{k) and tc{k— 1) to a linear function 

ti^\k)+aik)/k removes the leading l/k dependance and 
gives an even faster converging new sequence of points 
ti^\k). For the triangular lattice the result is shown in 
Fig.(|lO|). Whithout further extrapolation this confirms 
the above value of the critical point. 

Having located tc we now turn to a more detailed in- 
vestigation of the critical behaviour. Scaling theory pre- 
dicts this phase transition to be in the universality class 



of the three-dimensional xy-modelll^l. The Mott lobe is 
thus expected to close following a power law 

A{t) oc (tc - ty , (t,-t)«l . (18) 

where z the dynamical exponent is given by z-^ 1. From 
the scaling form of the dynamic susceptibilityH the power 
law divergence of the static structure factor can be de- 
rived 

Siq = 0, <) oc (tc " ty^ , (tc - t) « 1 (19) 
with 73 = (1 - ?7)iy (20) 

We used a Dlog-Pade analysis to check these scaling pre- 
dictions. For the square lattice this gives for the critical 
point tc ~ 0.0598 and the critical exponents « 0.69 
and 7s = 0.653. The different Fade approximants dif- 
fer quite a bit and the tc value clearly does not agree 
with the result of the 1/n extrapolation. The values for 
the critical exponent v have to he compared with the 
known value for the 3D xy-modelEj, j/ = 0.6693 ± 0.0010 
and 75 — 0.64.. ± 0.0022 , obtained by Borel summation 
of field theoretical results. Also the difference between 
the two results is quite small of the order of a few per 
cent this is still closer to the 3D Heisenberg model with 
v = 0.7054 zbO.OOlOy then to the xy-result. The conver- 
gence for the square lattice is not that good as mentioned 
before. When going to triangular lattice the accuracy of 
the analysis improves significantly. The Dlog-Pade ap- 
proximands find the critical point at tc ~ 0.03778 with 
exponent ly = 0.681 - 0.683 and 7^ = 0.651 - 0.652. The 
values for different approximants are presented in Tab. 
B. The tc value is in excellent agreement with the 1/n 
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FIG. 5. Square lattice: Phase Diagram obtained by Pade 
analysis. 



extrapolation. The critical exponent still remains higher 
then the field theoric value. Looking carefully at the 
values for tc and v as given by different approximants 
there is a tendency for both quantities to decrease slowly 
with increasing length of the series. In addition to the 
standard Dlog-Pade analysis we also performed a biased 
analysis with the critical point fixed at tc = 0.03778. 
This gave i' = 0.6805 - 0.6810. 

The critical exponent depends linearly on the critical 
point. Fixing v at its field theoretic value this would 
require tc = 0.03775 — 0.03776 which can clearly be ruled 
out by the l/n extrapolation. 

The discrepancy is probably due to subleading correc- 
tions to the asymptotic divergent term. The Dlog-Pade 
analysis ignores these contributions and assumes a pure 
power law divergence at the critical point. An improved 
analysis to resolve this subtle question however would 
require much longer series which won't be possible, be- 
cause the numerical problem grows exponentially with 
the order of the series. 



V. ID PHASE DIAGRAM 

We |iipw turn to the one-dimensional case. Scaling 
theoryllJ predicts the critical behavior of the system to 
be that of a Kosterlitz-Thouless transition for which the 
Mott lobes closes according to 

A(i)cx^expf--=^=) , (tKT-t)«l . (21) 

V V^KT — 1/ 

This highly nonanalytic behavior renders a direct series 
extrapolation complicated and the results of such a sim- 
pleminded approach are not conclusive. In order to de- 
termine the phase diagram we thus performed a Pade 
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FIG. 6. Triangular lattice: Phase Diagram obtained by 
Pade analysis. 



analysis for In^ A(t) cx (^kt — t)~^. This quantity has a 
simple pole at the critical point which can be captured 
by a rational function. 

The phase diagram obtained this way is shown in 
Fig.|ll| where we compare results from the series anal- 
ysis with numerical pdata of QMC simulations by, Ba. - 
trouni and ScalettaiB and a recent DMRG studyBEl. 
The agreement between the series and the DMRG data 
is excellent. Quite remarkable is the reentrance behavior 
found for t/U > 0.2. The series analysis thus confirms 
this phenomenon first observed in the DMRG calcula- 
tion. This type of behavior is a new feature not seen 
before in investigations of Bose systems. 

A simple intuitive way of understanding this phe- 
nomenon is the fact that Mott lobe is particle-hole asym- 
metric for the lattice problem. Starting from strong cou- 
pling {U ^ t) it is clear that the effect of the kinetic 
energy is to delocalize the particles. The delocalization 
decreases the average number of particles per site if the 
chemical potential is held fixed. On the other hand in 
the weak coupling limit (U ^ t) the bosons condense at 
the lower band edge so that for increasing bandwith (t) 
and fixed chemical potential the average number of parti- 
cles per site is decreasing. The nonmonotonic behavior of 
the density can be understood simply as a result of two 
limiting cases. Thus starting from the Mott phase the 
number of particles per site first decreases and then in- 
creases leading to a second Mott transition for a well de- 
fined range of the chemical potential, plot of n = const. 
lines in mean field ? 

The uncertainties in the precise location of the 
Kosterlitz-Thouless transition are still comparatively 
large. We use a Pade analysis of In^ A{t) cx (^kt — t)~^- 
This methods turned out to give excellent results. We es- 
timate the point for Kosterlitz-Thouless transition to be 
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FIG. 7. Square lattice: 1/k extrapolation of the critical 
point tc- The circles and squares are estimates from series 
for the gap and correlation function respectively. Note the 
excellent agreement between the tc estimates of these two 
independendant quantities. 



located at Ikt/U = 0.26±0.01 and hkt/U = 0.16±0.01. 



VI. CONCLUSION 

In conclusion, series expansion techniques were ap- 
plied to investigate the zero temperature properties of the 
Bose-Hubbard model. We determine the complete spec- 
trum of single-particle and single-hole excitations in the 
Mott phase. The phase diagram in one and two dimen- 
sions is obtained and the critical end points of the Mott 
insulator regions are determined. In two dimensions this 
is so far the only investigation of the complete phase dia- 
gram of this problem. In one dimensions the series shows 
almost perfect agreement with a recent DMRG study and 
provides a conclusive confirmation for reentrence behav- 
ior from the compressible to the insulating phase near 
the Kosterlitz-Thouless point. 

We achnowledge useful and interesting discussions 
on this problem with M. P. Gelfand, T. Giamarchi, 
T. Kiihner, A. J. Millis, R. Noack, A. v. Otterlo, 
R. R. P. Singh, U. SchoUwock, G. Schon, H. Schulz, and 
S. R. White 



APPENDIX A: EXPLICIT SERIES EXPANSION 
FOR SMALL CLUSTERS 



Perturbation theory for physical observables 
formulated as a coupled cluster expansion, see e.g 
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FIG. 8. Square lattice: 1/n extrapolation of the critical 
point sequence ti^' . Note the scale of the vertical axis! 



allows to write the ground state energy for the Bose- 
Hubbard model on a particular lattice C in the following 
way 



EMottiC)/N^E^2jN + y^Lig)xWEig) ■ (Al) 



Here g are the connected clusters or graphs that can be 
embedded into the lattice £, the so called lattice constant 
L(g) is the number of embeddings per site of g in £ and 
Wsig), the cluster weight, is the contribution of g to the 
expansion of the energy E. Each cluster represents itself 
a small lattice and the weights can be expressed by the 
recursion relation 



WEig) = Emottig) - ^ WE{g') 



(A2) 



where i^Mottls) is the ground state energy of the par- 
ticular cluster g under consideration and the g' are all 
its subclusters. An important property of the weights is 
that 



WE{g)^0{x^) 



(A3) 



for a cluster with n bonds (x — t/U). The procedure 
to calculate the series for the ground state energy up to 
a certain order 0{x'^) is now straightforward: first one 
has to find all cluster with up to n bonds and determine 
their lattice constants L{g), in the next step for each 
cluster its ground state energy i?Mott(5) has to be deter- 
mined perturbatively from which the weights follow by 
eqn. ( |A2| ) . Summing up all the contributions, eqn.(Al) 
gives the desired result. 

As an example we briefly discuss the expansion for the 
Bose-Hubbard chain. The clusters that contribute are 
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FIG. 9. Triangular lattice: 1/fc extrapolation of the critical 
point tc- The circles and squares are estimates from series 
for the gap and correlation function respectively. Note the 
excellent agreement between the tc estimates of these two 
independendant quantities. 



particular simple: finite open chains with n bonds and 
n+1 sites all having lattice constant L{g) = 1. The first 
three cluster that contribute have the following ground 
state energies: 



E{9) 
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Nest step is the subcluster subtraction: the one bond 
cluster can be embedded two (three) times into two 
(three) bond cluster and the two bond cluster also two 
times into the three bond cluster. This gives for the 
weights: 
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FIG. 10. Triangular lattice: 1/n extrapolation of the criti- 
cal point sequence ti^\ Note the scale of the vertical axis! 



the ground state energy: 

272 

Emou {C = 1D)/N = -4 x2 -I- 4 x^ -I- — x^ 

y 

We now turn to the construction of the effective hamil- 
tonian for the excited states. As an example we will 
discuss the single particle states, eqn.(||). The effective 
hamiltonians obtained by degenerate perturbatition the- 
ory for the clusters with up to three bonds are given be- 
low. Note that the energies are given with respect to the 
ground state energy, because the coupled cluster theorem 
applies only to the difference Hcs — EqU. 
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Summing up all these different term gives the series for 
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FIG. 11. Phase Diagram of the Bose-Hubbawl chain: Nu- 
merical data from DMRG r^Ught circles, ref.(y) and QMC 
studies (dark squares, ref. (□) compared to the results of a 
Pade analysis. 
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The cluster hamiltonians are not hermitian in all orders 
of the expansions. While this looks strange at first sight it 
is just a consequence of the finite cluster size which lacks 
the translational symmetry of the full lattice. Therefore 
the sites of a finite cluster arc no longer all equal and so 
the sequence of intermediate virtual states that connects 
two different sites i and j depends on the actual order, 
i.e. in perturbation theory going from site i to site j 
is not equivalent to going from j to i. An example is 
the following sequence of states that contributes in third 
order to H{1, 2) the hamiltonian of the two bond cluster: 



|2, 1, 1 > 



2,2,0> 



1, 3, > 



1,2,1> 



In the contribution to i?(2, 1) the same intermediate 
states will contribute but in reverse order and this gives 
rise to the non hermitian cluster hamiltonian. There is 
no equivalent sequence to the one shown above, because 
the sites 1 and 2 are not equivalent. When the effective 
hamiltonian for the full lattice will calculated later on 
all embeddings of the cluster into the lattice have to be 
taken into account. This restores the the translational 
symmetry and the resulting Hamiltonian will be hermi- 
tian again. 

A new complication arises when performing the sub- 
cluster subtraction, because matrix elements have to be 
subtracted and it is therefore necessary to keep track of 
exactly which sites of a cluster are covered by the sites 
in a particular subcluster. The subtraction procedure is 
shown in the following graph. 
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1 

•- 



2 
-•- 



1 

X— 

1 

(1) 

1 



2 

2 

(1) 



3 

3 

(1) 

3 
O 



(1) (2) 



1 
O 



2 

•- 



(1) (2) 



The subtraction procedure gives the cluster weights of 
the hamiltonian and as it is the case for the ground state 
energy the weights for a cluster with n-bonds are at least 
of order 0(a;"). 



1 



1 

•- 



Wh{1,1) = 1 



2 



Wif(l,2) 
Wh (2,1) 
Wh (2,2) 

2 3 



Oa; 
2x 
2x 
Ox 



5/2 x^ 



Ox^ - 
Ox^ - 



3/2 x^ 
3/2 x^ 



5/2 x^ 



Wh{1,2) 


= - 


\-Qx + 


Ox^ - 


9a;3 


Wh{1,3) 


= - 


fOa; - 


2x^ + 


Qx^ 


Wh {2,1) 


= - 


\-0x + 


Qx^ + 


3/2 a;3 


Wh{2,3) 


= - 


\-0x + 


Ox^ + 


3/2a;3 


Wh{3,1) 


= - 


fOa; - 


2x^ + 


Qx^ 


Wif(3,2) 


= ' 


f Oa; + 


Qx^ - 


9a;3 



Wh{IA) = +0a; 
Wif(4,l) = + 0a; 



Qx^ 
Qx^ 



6a;3 
Qx^ 



The last steps requires to find all possible embeddings 
of the cluster in the chain. Due to the translational in- 
variance of the problem, it is sufficient to search only 
those possibilities where that connect the origin to other 
sites. For the first clusters this is demonstrated in the 
following: 
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Up to third order in the expansion parameter the ef- 
fective hamiltonian for the single particle states in the 
Bose-Hubbard chain reads: 



Hcsi~3) 


= - 


fOa; 


+ 0x^ 


+ 6x^ 


i^eff("2) 


= - 


fOa; 


~2x^ 


+ 0x^ 


Hesi-1) 


= - 


f 2a; 


+ 0x^ 


-9x^ 


iJeff(O) 


= f - 


fOa; 


+ 5x^ 


+ 6x^ 


i?eff(l) 


= - 


f 2a; 


+ 0a;2 


-9x^ 


i?cff(2) 


= - 


fOa; 


~2x^ 


+ 0x^ 


i?cff(3) 


= - 


fOa; 


+ 0x^ 


+ 6x^ 



Note that the effective hamiltonian is hermitian. 



APPENDIX B: HASH FUNCTION FOR LATTICE 
BOSONS 

Perturbation theory requires to work in the full hilbert 
space of the finite clusters. It is therefore necessary to 
have an efficient method to represent the basis and the 
hamiltonian. In particular there is a need for a Hash 
function that labels each basis state in a unique way. 
The hamiltonian (|l]) conserves the total particle number 
N. For a cluster with L sites and fixed iV the hilbert 
space has a total of 



D{N,L) 



1)! 



m {L-iy. 



(Bl) 



states. The natural basis is given by the set of occupation 
numbers {ui} of all lattice sites | ni , n2 , na , ... , > 
. We now describe how to construct the full basis. The 
first state is the one where all particles are sitting on site 
1. Then follow those states were the first two sites are 
occupied. In the next stepthe last occupied site is i = 3. 
From those possibilities we first get all those with just 
one boson on site 3. There are then three bosons left 
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to be distributed on i = 1 and i = 2 and we repeat the 
procedure, first all three on site 1 and so on. 
Define 



:=^n, , 00 -.= 



(B2) 



Which is related to the occupation number by 

n, = e,, - 6,^1 (B3) 

Hash function 

L 

Hf{{m})=J2 [D{d,,i)-D{9,_,,i)] (B4) 
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TABLE I. Scries for the single particle gap A annd the 
equal time structure factor 5q=o of the Bose-Hubbard model 



TABLE 11. Scries for the single particle gap A armd the 
equal time structure factor 5q=o of the Bose-Hubbard model 



on 


the square lattice. 


on 


the triangular lattice. 




gap series 




gap series 


n 




n 


On 





1 





1 


1 


-12 


1 


-18 


2 


-22 


2 


-81 


3 


-264 


3 


-819 


4 


-15659 / 10 


4 


-54891 / 4 


5 


-656984 / 25 


5 


-11459377/50 


6 


-513092341 / 2250 


6 


-6764830501 / 1500 


7 


-13396365654 / 3375 


7 


-3957593443549 / 45000 


8 


-2194497431888101 / 56700000 


8 


-23724030434424597 / 12600000 


9 


-523244582353596437 / 744187500 


9 


-79054659543137812691 / 1984500000 


10 


-4749112579154967367231 / 625117500000 


10 


-49507676563116513700717 / 55566000000 


11 


-6676218845916748474723399 / 49228003125000 


11 


-19752788544107.0312336222045518956 


12 


-5669114326328304841982042447 / 3508972062750000 


12 


-454652221307484.008486685913608139 


13 


-3473317126780784521271398100 / 126449317253259 


13 


-10398395856680122.5582305335816799 




structure factor series 




structure factor series 


n 


On 


n 


On 





3 





3 


1 


32 


1 


48 


2 


432 


2 


1080 


3 


6656 


3 


25344 


4 


99632 


4 


613016 


5 


14154496 / 9 


5 


15108128 


6 


663550400 / 27 


6 


3391856144/9 


7 


31905307840 / 81 


7 


28444402112 /3 


8 


7618958766796 / 1215 


8 


97228564590772 /405 


9 


12934227681606432 / 127575 


9 


86595555599744452 / 14175 


10 


14570617373829713351 / 8930250 


10 


2787620342817465447617 / 17860500 


11 


19905912307529372064253 / 750141000 


11 


82527451969123616224435919 / 20628877500 


12 


17231660529006598072716626038 / 40068302174901 


12 


102821325219551430.224270929006395 


13 


47514032474492554578981737799.4028 / 6764778289269 


13 


2648736908451106507.61278132531748 
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TABLE III. Series for the single particle gap A aniid the 
equal time structure factor S'q=o of the Bose-Hubbard chain. 



gap series 


n 


an 





1 


1 


-6 


2 


5 


3 


6 


4 


287/20 


5 


17463 / 150 


6 


-1806729 / 3000 


7 


73674531 / 22500 


8 


-297690613629 / 16200000 


9 


14666046468323 / 121500000 


10 


-6295148943458549 / 7290000000 


11 


114441271150219589 / 18225000000 


12 


-422231271662550684871 / 9185400000000 


13 


1473292023890353319230511 / 4340101500000000 


structure" factor sc^rics 


11 


tin 





3 


1 


16 


2 


72 


3 


320 


4 


4120/3 


5 


48544 / 9 


6 


596992 / 27 


7 


2396512/27 


8 


27653840 / 81 


9 


1329678 


10 


22598877209 / 4375 


11 


1275342277201 / 65610 


12 


4055776421430107/55112400 


13 


203597243119484303 / 723350250 



TABLE IV. D-log Pade approximants for the critical point 
tc and the critical exponents v and 7s of the Bose-Hubbard 
model on the square lattice. 





gap 


structure factor 


[L/M] 




tc 








tc 




Is 


[4/4] 





05980859 





69137519 





05976303 





65480784 


[4/5] 





05978709 





68924487 





05976978 





65534864 


[5/4] 





05980549 





69112808 





05971761 





65393342 


[4/6] 





05980887 





69140624 





05976747 





65515773 


[5/5] 





05980805 





69133305 





05976666 





65508534 


[6/4] 





05980880 





69140183 





05976814 





65522423 


[4/7] 





05980331 





69097429 





05975849 





65433135 


[5/6] 





05984310 





69384313 





05978395 





65596948 


[6/5] 





05980488 





69108539 





05976313 





65477892 


[7/4] 





05980859 





69138480 





05976778 





65519366 


[4/8] 





05978512 





68975552 





05974697 





65294432 


[5/7] 





05987138 





69556354 





05973294 





65050768 


[6/6] 





05987084 





69552706 





05972443 





64864666 


[7/5] 





05978277 





68960137 





05975399 





65386465 


[8/4] 





06016134 





72645361 





05974142 





65212846 



TABLE V. D-log Pade approximants for the critical point 
tc and the critical exponents v and 7s of the Bose-Hubbard 
model on the triangular lattice. 





gap 


structure factor 


[L/M] 




tc 




V 




tc 




Is 


[4/4] 





03781415 





68753874 





03768377 





60598248 


[4/5] 





03782464 





68812741 





03778221 





65202471 


[5/4] 





03777127 





67855607 





03778121 





65179657 


[4/6] 





03780279 





68619855 





03777688 





65066891 


[5/5] 





03780194 





68605142 





03777941 





65137574 


[6/4] 





03780562 





68666551 





03777913 





65130264 


[4/7] 





03779156 





68397787 





03777917 





65130976 


[5/6] 





03776841 





67611321 





03777930 





65134670 


[6/5] 





03779251 





68419813 





03777930 





65134598 


[7/4] 





03778693 





68276629 





03777931 





65134883 


[4/8] 





03778359 





68180328 





03777864 





65116474 


[5/7] 





03778045 





68074745 





03777854 





65113450 


[6/6] 





03778034 





68070944 





03777943 





65137873 


[7/5] 





03778358 





68180026 





03777915 





65130775 


[8/4] 





03778162 





68114199 





03777683 





65058408 
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